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Abstract 

|/-\ . We propose a conceptually and computationally simple method to evaluate the neutrinos emitted by 

£N| supernova remnants using the observed 7-ray spectrum. The proposed method does not require any prelim- 

inary parametrization of the gamma ray flux; the gamma ray data can be used as an input. In this way, 
(-H ' we are able to propagate easily the observational errors and to understand how well the neutrino flux and 

O |. the signal in neutrino telescopes can be constrained by 7-ray data. We discuss the various possible sources 

JL ' of theoretical and systematical uncertainties (e.g., neutrino oscillation parameters, hadronic modeling, etc.), 

j_( . obtaining an estimate of the accuracy of our calculation. Furthermore, we apply our approach to the su- 

(y>j ' pernova remnant RX JI7I3. 7-3946, showing that neutrino emission is very-well constrained by the H.E.S.S. 

7-ray data: indeed, the accuracy of our prediction is limited by theoretical uncertainties. Neutrinos from 
RX J1713. 7-3946 can be detected with an exposure of the order km 2 x year, provided that the detection 
threshold in future neutrino telescopes will be equal to about 1 TeV. 
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1 Introduction 

Under-water and under-ice neutrino telescopes are instruments aiming to discoveries. They could reveal 
an effective acceleration of cosmic rays (CR) in galactic sources (such a supernova remnants pQ and 
micro-quasars [5]) and/or in extragalactic sources (such as AGN and gamma ray bursts [H])- However, at 
present it is possible to obtain reliable expectations only for few of these sources, such as the supernova 
remnants (SNR) discussed in this paper. 

The idea that CR could originate in supernovae has been put forward already in 1934 but the first 
quantitative formulation of the conjecture that the young SNRs refurnish the Milky Way of cosmic 
rays, compensating the energy losses, is due to Ginzburg & Syrovatskii pQ. In fact, the turbulent 
gas of SNRs is a large reservoir of kinetic energy and this environment can support diffusive shock 
waves acceleration |3]. The theory of CR acceleration in SNRs is still in evolution, but the generic 
expectations are stable. The CR flux in SNRs is expected to have a power law spectrum with spectral 
index T = 2.0 — 2.4 at low energies with a cutoff at an energy E c which depends on the details of the 
acceleration mechanism and on the age of the system. In specific implementations the cutoff energy 
E c can be as large as several PeV and, thus, is consistent with the "knee" in the CR spectrum at 
E ~ 3 x 10 15 eV, which is believed to mark the transition from galactic to extra-galactic origin of 
CR 0. 

In recent times, great progress has been made in the observation of SNRs. In particular, the High 
Energy Stereoscopic System (H.E.S.S.) [6] has determined quite precisely the gamma ray spectra of few 
SNRs showing that they extend above 10 TeV. In the context of Ginzburg & Syrovatskii hypothesis, 
it is natural to postulate that the observed gammas are produced by the decay of it (and rj) resulting 
from the collision of accelerated hadrons with the ambient medium. New and crucial observations are 
being collected and the hadronic origin seems to be favored for certain SNRs, such as Vela Jr [7] and RX 
J1713. 7-3946 [El[9]|j It is not yet possible, however, to exclude that (part of) the observed radiation is 
produced by electromagnetic processes. The definitive proof that SNRs effectively accelerate CR could 
be obtained by the observation of high energy neutrinos in neutrino telescopes presently in operation, 
in construction or in project |llj . 

As well known, there is a strict connection between photon and neutrino fluxes produced by hadronic 
processes in transparent sources (see, e.g., [12]), which results from the fact that the same amount of 
energy is roughly given to ir°, 7r + and tt~ in hadronic collisions. On this basis, one can estimate the 
neutrino fluxes expected from a sources with known 7-ray spectra, trying to identify detectable sources 
and/or to optimize the detection strategies. The gamma-neutrino connection has been described in 
various recent papers [HI HU HU HU [THJ, HJ EUJ, [21] at a different level of accuracy, relying on different 
assumptions on the primary cosmic ray spectrum and/or hadronic interaction model. 

In this particular moment, when high energy gamma ray astronomy is flourishing and the neutrino 
telescopes are finally becoming a reality, it is clearly important to have solid and transparent predictions. 
We continue, thus, the work started in [16\ [2~T] 122] , proposing a method to calculate neutrinos fluxes 
which is, at the same time, simple, accurate and model-independent. Our results are in essence a 
straightforward applications of standard techniques [23], but we believe that that they will be useful 
since they improve the existing calculations in various respects. In particular, we provide simple analytic 
expressions for the neutrino fluxes which have a general validity and can be applied directly to gamma 
ray data since they do not require any parametrization of the photon spectrum. This allows us to 
propagate easily the observational errors in the gamma ray flux and, thus, to understand how well 



1 See also [10] for a recent analysis leading to a different conclusion. 



the neutrino flux and the signal in neutrino telescopes can be constrained by 7-ray data. We also 
discuss the various possible sources of theoretical and systematical uncertainties, obtaining an estimate 
of the accuracy of our method. This kind of analysis is relevant in the present situation, since the 
number of events expected in neutrino telescopes from SNRs is quite low (see, e.g., [16]) and even 
small (downward) revisions of the expected signals may be important and/or require different detection 
strategies. It is thus important to understand the relevance of different assumptions in the calculations 
and the origin of (apparently) contrasting results appeared in the literature. 

The plan of the paper is as follow. In Sect. [2] we review the method presented in our previous 
paper [16]. The relations obtained in this work - Eqs. ([8]) - are physically equivalent to those presented 
in [16]. However, they are more compact and more convenient for numerical computations since we 
have been able to recast them in such a way that they require only one numerical integration. In 
Sect. [31 we discuss the effect of neutrinos oscillations (see [24] for a review) and we discuss the relevance 
of uncertainty in neutrino mixing parameters for the predicted neutrino flux. In Sect. H] we present 
our main results, i.e., Eqs. (|30p and (|32|) which relate the (oscillated) muon neutrino and antineutrino 
fluxes to the 7-ray flux. As it is explained, the only necessary information to predict neutrino fluxes are 
the relative production rates of the various mesons in hadronic processes, which are robust predictions 
of hadronic interaction models. In this paper, we adopt the results from Pythia [25], estimated by 
using the parametrization of hadronic cross sections presented in [26] . A comparison with SYBILL |27| 
and DPMJET-III [28] is performed (when possible) by using the parametrization and tabulations of 
hadronic cross sections presented by [29] and [30]. In Sect. [5] and [6] we propose a procedure to predict 
neutrino fluxes and event rate in z^-telescopes directly from 7-ray observational data, and we apply it to 
RX J1713. 7-3946 which is presently the best studied SNR. In Sect. [7] we summarize our main results. 

2 The photon-neutrino relation 

The interactions of cosmic ray protons with a hydrogen ambient cloud result in the production of mesons 
which subsequently decay producing gamma ray and neutrinos. Both gamma rays and neutrinos depend 
linearly on the flux of the primary cosmic ray. We, thus, expect that a linear relation also exists between 
photon and neutrino fluxes, which can formally be expressed as: 






as will be precised just below in this Section. 

In order to determine the kernels K V \E V ,E^\, we only need to know the relative number of pions, 
kaons and r\ produced by cosmic rays at any given energy, as it is explained in [16] and further discussed 
in the following. Let us indicate with Ri[E] the number of i— particles produced per unit time and unit 
energy in the cloud. Here, for simplicity, we assume that the momentum distributions of the various 
particles are approximately isotropic, so that the differential flux produced in a detector at a distance 



D (if the i— particle is stable) is simply given byo 

Photons are mainly produced by 7r° and 77, according to: 

2 Z" 00 dE 

where 6 r)7 = 0.394 is the 77 — > 77 branching ratio and the factor 2 takes into account that two photons 
are produced for each 7r° and 77. Neutrinos are, instead, mainly produced by charged pions and charged 
kaons. We can formally write (neglecting neutrino oscillations): 

1 /-oo at? 

® AK] = 4^^2j E ^RilEMK/E] (3) 

where i = tt + , ir~ , K + , K~ and v = u e ,u e , v^, V^. The quantity Wi u [x]dx, with x = E u /E, represents 
the spectrum of neutrinos v produced in the decay chain the i— meson. 

If we assume that the ratios between the production rates of the various mesons are approximately 
constant (see next section), the above relations can be combined in order to obtain the neutrino fluxes 
as a function of the photon flux. Rel. (|2|) can be, in fact, inverted obtaining: 

. . 4irD 2 E d$JE] 

RAE] = -TTKj n l^ (4) 

where f v = (R^/R^o). This expression can then be used in Eq. ([3]) to obtain: 

POO AT? 

ME A = -E^K V [E v /Ey] $ 7 [£ 7 ] (5) 

JE V &j 
where: 



K \ x] - I y f ^iM (q) 

u[X] ~ 2{l + b rn f ri )^ U d\nx W 

and 

f l = R l /R n o wlthi = Tr+,TT-,K + ,K-, V . (7) 

In explicit terms, we can write: 

Kv, [ x ] = T—77 (/tt+ g-xve [x] + I'k+ 9Ku e [x] ) 
Kv e [x] = —-JJ (/*- 9iru m [x] + f' K - 9Ku e [x]) 

K Vfk [x] = (f w + K [x] + U- ftn,„ [x] + f' K+ h K [x] + f' K - g Kv>l [x]) 

K V» [x] = ——77 (L+ ftri/ p [x] + L- K [x] + f' K + gKv^ [x] + f K - h K [x]) (8) 



If this assumption is removed, one has to replace here and in the following: 

l -RAE] - dR ^ n] 



4tt l ' dQ 

where dRi/dD. is the rate of i— particles produced per unit energy and unit solid angle, n is the unit vector in the direction 
connecting the SNR to the detector and we have taken into account that the produced particles are almost collinear. 



where we have defined 

f'r, = b vifv /x± = b Kufx± (9) 

and bKv = 0.634 is the branching ratio for K^ semi-leptonic decay. The functions hi[x] account for 
neutrinos produced in 7r + — > fj, + + v^ and K + — > fj, + + v^ (and charge conjugated processes) and are 
given by 

h l [x] = ^5[x-(l-r l )] (10) 

where 7*j = (ttt-^/ttij) 2 with i = it, K. The functions <fo,[x] account, instead, for neutrinos produced 
by muons. They also encode the information on the energy distribution and polarization of muons 
produced by pions and kaons decay and can be expressed as: 
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(12) 



(13) 

Some details of the derivations and some checks of these formulae, in the limit of "isospin invariance" , 
are given in the appendix A of |21j . 

We note that the same kind of approach can be used to calculate the flux of electrons and positrons 
•produced by hadronic processes in the cloud, which may be relevant to calculate radio synchrotron 
emission. By considering that the spectral distribution of electrons produced in fi — ► e + V e -\- v„ is 
identical to that of of i^, we immediately obtain: 



**±[E\= ^K e± [E/E,]^[E,} (14) 



where 



K e~ [ x ] = J— 77 {fir- 9™» M + f' K - 9Ku„ [x]) (15) 

K e+ [x] = — -JJ {L+ 9nu„ [x] + f' K + 9Kv„ [x]) (16) 

3 The effect of neutrino oscillations 

Neutrino telescopes are sensitive to muon neutrino and muon antineutrino fluxes at earth that differ 
from the fluxes produced in the supernova remnants due to the effect of neutrino oscillations. Denoting 
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Figure 1: Suppression of cosmic v^ (or v^) flux as a function of the electron /muon neutrino (or antineutrino) 
flux ratio. The dashed (continuous) lines enclose the la (2a) region consistent with the measurements. The 
dotted line identify the most probable value for the suppression factor. The grey region, instead, is forbidden. 



Puu$L + Peu<£>v e , where due 



with <&° (£ = e, fi, t) the flux in absence of oscillations, we have <J> t 
to cosmic distances, the oscillation probabilities are given by the formula Pee' = J2i=i 2 3 Wu\We'i\, 
appropriate for averaged oscillations [31] . This formula is symmetric in the exchange £ «-> £' and is valid 
for neutrinos and antineutrinos. 

In order to implement oscillations in our approach, we replace the kernel in Eq. © with the the 
"oscillated" kernel: 
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(17) 
(18) 



which is given by: 

K™ [x] = P^ K V)1 [x] + P eil K Ve [x] . 

This can be expressed as the sum of several contributions, each arising from a different mechanism for 
neutrino production: 
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- 1 *■ 1/ 



A\ h n [x] + A 2 gnu,, [x] + A 3 g nV£ [x] + 
+ -Bi h K [x] + B 2 g K vn. [x] + B 3 g Kve [x] 



(19) 



with the coefficients Ai and Bi given by: 
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(20) 



A similar expression holds for the muon antineutrino flux. Namely, the kernel i^ sc [x] is obtained from 
Eq. (fT9|) with the replacements /„.+ — > /„.- and fx+ — > fx- (and vice versa) in Rels. (f20|) 

The oscillation probabilities P^ and P eil can be evaluated by assuming: 0\ 2 = 35° ±1°, 823 = 42° ±4° 
and #13 = 5° ± 4° , close to the range found in a recent global analysis of the world data [32] . We obtain 
the central values: 
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0.36 



P 



e/j 



0.26 



(21) 



that will be used in the following calculations. The uncertainties in P^ and P efjL are at the level of 
about ~ 10% and are almost completely anti-correlated, being mostly due to the spread of #23 around 



maximal mixing. It is, thus, not immediate to understand how they propagate to the final muon 
neutrino flux. Writing: 

% = — ^-(Puu + P e „V) where i e M (22) 

we see the effect of oscillations depends on the neutrino flavor ratio at production ?/>. We also understand 
that we expect a partial cancellation between the contributions of P w and P efM to the total error budget. 

In Fig. [H we show the l.h.s. of the above equation as a function of ip. The dashed (continuous) lines 
in the figure enclose the la (2a) region consistent with the measurements. The dotted line identifies 
the most probable value for the muon neutrino flux suppression factor. Several remarks on this figure 
are in order. 

(i) When we consider the values ip ~ 0.5 characteristic of neutrinos from pion decay, we find that 
the most probable value is $j, M /<3?* ~ 0.33. This is a well known result, first derived in [33] and often 
used in the literature (see, e.g., [HI EI]). For the sake of precision, we note that ip is the electron/muon 
neutrino (or antineutrino) flux ratio at a fixed energy which may differ from 0.5 even for pion decay, 
due the different energy distribution of the three neutrinos produced by pions. If we assume that the 
neutrino spectra in the SNR is described by a power law with spectral index a = 2 (a = 3), we obtain 
ip = 0.54 (ip = 0.61) o In our approach, however, we do not need to decide a priori the neutrino flavor 
ratio ip or the value of <&j/ /<£* ot : the energy distributions of the produced neutrinos and the oscillation 
effects are automatically implemented by Eq. (fl~9|) . 

(ii) We note that the possible suppression of the muon neutrino flux is bounded from below, as 
emphasized by the forbidden (grey) region in the figure. In the extreme case ip = we have $j, /$* ot > 
1/3 that can be understood by considering that the (averaged) survival probability P„^ cannot be 
smaller than 1/3; if we assume ip = 1 we obtain, instead 3V M /<3?* ot = (P^ + P e ^)/2 = (l-P flT )/2 > 1/4, 
which follows from the fact that the (averaged) oscillation probability P^ r cannot be larger than 1/2. 
The oscillation effect that is realized in Nature happens to be very close to the maximum possible effect, 
namely, to be very close to the forbidden region in Fig. [TJ 

(in) The uncertainty due to the imprecise knowledge of the oscillation parameters is at the level of 
2% and diminishes with increasing ip, as a result of the partial cancellation of the much larger (anti- 
correlated) contribution of P w and P efl to the total error budget. It should be noted that uncertainty 
works mostly in the direction to increase the expected muon neutrino (or antineutrino) flux. For 
ip ~ 0.5, we find, in fact, the region 

|^ = (0.33 - 0.35) at 2a (23) 

which is extremely asymmetric with respect to the most probable value ~ 0.33. 

4 The meson production rates (the fa factors) 

The last information that we need to predict the neutrino fluxes are the factors fi which give the 
production rates of the various mesons rescaled to that of the neutral pions (see Eq. ([7])) and can be 
calculated from hadronic interaction models. 



3 The neutrino flavor ratio ip for K decay can largely depart from 0.5. Assuming, again, that the neutrino spectra is 
described by a power law with spectral index a — 2 (a = 3), we obtain ip — 0.32 (ip = 0.19). In a more realistic case, in 
which we choose / w + = f n - = 1 and f K + = f K - = 0.1, we obtain ip = 0.52 (ip = 0.56). 



The rate of production of i— mesons at an energy E is given by: 



Ri[E] = N 



dE v 
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(24) 



where N is the total amount of target hydrogen in the cloud, Jp^p] is the CR energy spectrum (averaged 
over the SNR volume, see Eq. (7) in [22]), and we employed the usual definition of the adimensional 
distribution function Fi[x,E p ]: 

duj _ cr[E p ] 

dE " 
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(25) 



where a is the total inelastic p-p cross section, E p is the proton energy, dai/dE is the inclusive cross 
section for i— particles production and E represents the energy of the produced particle. Let us assume 
that the cosmic ray spectrum is roughly described by a power law in energy. We can write J[£ , p ]cr[£'p] oc 
E~ a , and using the "quasi-scaling" approximatiorQ for hadronic cross sections, i.e., taking 



where E p is a fixed reference values for the proton energy, we finally obtain: 



fi 



Ri 



Zj[a] 
R^o Z n o [a] 



where: 



Zi[a] = f dxx a ~ 1 F i [x] 
Jo 



(26) 



(27) 



(28) 



The above relation shows that the ratios between the number of produced particles are essentially 
independent on energy and are determined by the (a — 1)— momenta of the adimensional distribution 
function Fj. 

Presently there are various well developed codes - such as Pythia [25], SYBILL [27], QGSJET [35] . 
DPMJET-III [28] - which implement different hadronic interaction models to study particle production 
in nucleon-nucleon or nucleon- nucleus collisions. It would be interesting to make a systematic com- 
parison between the various codes, but this is beyond the scope of this work. In Ref. [26], suitable 
parametrizations of the pion and kaon differential cross section in pp collisions from the Pythia code 
are given. Moreover, parametrizations of 7r° and r\ productions in pp collisions, according SYBILL and 
QGSJET code, are presented in Ref. [29] • In the left panel of Fig. [21 we use these parametrizations to 
calculate the relevant fi factors for a = 2 — 4. We see that the quantities fi are marginally dependent 
on the assumed spectral index and we take the values obtained for a = 2: 



/tf-t 



1.08 
0.13 



f K - 



0.79 
= 0.09 



/„ = 0.48 



(29) 



as reference values in our calculations. We remark that, using these values, we improve on our previous 
calculation [16] where we assumed "isospin invariance" (i.e., f w + = f w - = 1) and we overestimated the 
rate of K^ production (/#-+ = fx- = 0-2). 



4 The "quasi-scaling" approximation is accurate at the few per cent level in the energy range of interest, see, e.g. 
We also remind that the total cross-section <t[B p ] has a very weak dependence on the proton energy, see, e.g., [29]). 



[22] 




Figure 2: Left panel: The factors / ff + (dashed blue), f„- (dotted blue), (/„.+ + / 7r -)/2 (so/id blue), f n (red 
line), fx+ (dashed black), f K - (dotted black), (fx+ + fi<-)/2 (solid black) calculated as a function of a according 
to Eq. p?7p . Right Panel: The kernels k Ufi [x] (solid) and k$ [x] (dashed). 



In order to estimate uncertainties in our method, we should compare the f\ factors in Eq. ()29[) with 
those obtained by using different hadronic interaction models. For assumed spectral index a = 2, the 
results of [30], which describes the various possible contribution to 7-production in p-ISM collisions, 
can be used to estimate the fi factors predicted by DPMJET-III. In this case, in fact, the factor 
Zi = (Ei)/E p coincides with the average fraction of the parent particle energy carried by particles of 
the type i. This quantity can be related to the energy fraction emitted in photons, by taking into 
account the branching ratios and the kinematic of the decay processes. By using this approach, we 
estimate (when possible) the values presented in Tab. Ql where we also show the predictions obtained 
from Pythia and SYBILL parametrizations. We remark the Pythia and SYBILL results are obtained by 
assuming pp collision, while the DMPJET-III results refer to the collision of protons with ISM, which 
is assumed to be composed by 90% protons, 10% helium nuclei, 0.02% carbon and 0.04% oxygen. One 
sees from Tab. [I] that there is a reasonable agreement between the prediction of the different codes. In 
particular, the fi factors are practically unchanged in the various cases. This means that the neutrino- 
photon ratio is a rather solid prediction, which should not suffer from large uncertainties in hadronic 
models and/or in the modelization of the CR flux and the ambient medium in the SNR. 

Finally, we discuss the other assumptions in our approach. Our calculation neglects the contribution 
of Kg and K® (and more rare production channels) to photon and neutrino production. According to 
Tab. [H these processes accounts for several percents increase of the photon and neutrino production 
rates with compensating effects in the neutrino-photon ratio, which is affected at the few percents level. 
Moreover, we assumed that the ratio between the production rates of the various mesons does not 
depend on energy, nor on the CR spectral shape. This assumption is motivated by the "quasi-scaling" 
behavior of hadronic cross sections and by the results presented in Fig. [2] and it is valid with few 
percents accuracy. In conclusion, by taking into account the approximations implicit in our method 
(e.g., constant fi factors, neglected production channels, etc.), the uncertainties in hadronic modeling 
and the uncertainty in neutrino oscillation parameters, we can safely estimate that the neutrino fluxes 
predicted by our approach are accurate at the level of about 20%|£| 



5 As recalled in the Introduction, a possible limit of our (and other) calculations is the possible presence of a relevant 
leptonic contamination (component) of the gamma radiation, that requires efforts in theoretical modeling of the sources 
and multi-wavelength observations to be precisely assessed. Such a possibility implies that, conservatively, we should 



Table 1: Prediction from different interaction models. See text for details. 
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Huang et al. [30] 
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Zi 
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0.16 
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Kelner et al. [29 

pp - SYBILL 



0.13 
1 



0.062 
0.48 



pp - SYBILL fi 1 0.48 

' Estimated by assuming that "direct 7 production" in [30] is due to n decays 



5 The neutrino flux 

By using the /, factors given in Eq. (|29p and the oscillation probabilities given in Eq. (|2ip . one can 
calculate the photon- neutrino kernels X° sc [x] and i^ sc [x] according to Eqs. (|19I20|) . Considering the 
explicit form of the the various contributions (see Eqs. (llO)lip ). one obtains the following simple analytic 
expression: 



$^[E] =0.380 ^[E/(l-r n )] + 0.0130 $ 7 [£/(l - r^)] + / — fe^ [x] $ 7 [£7/a 

Jo x 



(30) 



where the first two terms describe neutrinos produced in pions (first term) and kaons (second) decays. 
The kernel k u [x] , which takes into account neutrinos produced by muon decay, is shown in the right 
panel of Fig. [2] and it is given by: 



K, [ X ] 



x 2 (15.34- 28.93 x) 



0.0165 + 0.1193z + 3.747x 2 - 3.981x 3 



= (1 - xf (-0.6698 + 6.588x) 
A similar relation holds for muon antineutrino flux: 



x < r K = 0.0458 

Tk < X < 7V 

x > r n = 0.573 



(31) 



/■i dx 
^[E] = 0.278 $ 7 [£/(l - r n )} + 0.0090^^/(1 - r K )] + / — k Vfi [x] ^[E/x] (32) 

JO x 



where 



kp^x] = x 2 (18.48- 25.33 x) 

= 0.0251 + 0.0826x + 3.697x 2 - 3.548x 3 
= (1 - xf (0.0351 + 5.864x) 



x <r K = 0.0458 
rK < x < 7v 

X > Tir = 0.573 



(33) 



Eqs. pop and (|32p have a general validity and do not require any specific parametrization of the photon 
spectrum. They will be used in the following to derive the neutrino flux from the SNR RXJ1713. 7-3946 
directly from the observational 7— ray data. It is, anyhow, interesting to discuss the application of these 
relations for the functional forms most commonly used to parametrize the photon spectrum. 



speak of an upper bound on the neutrino signal within the present information. 
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5.1 Using parametrized fluxes 

If the photon spectrum is approximated with a power law <J> 7 oc E , one immediately sees that: 

$„[£] = Z V [V] ■ %[E] (34) 

where: 

Z [T] = 0.380 xl + 0.0130 x T K + / dx k Uu [x] x r ~ 1 

Jo 

Z Vti [T] = 0.278 x\ + 0.0090 x v K + J dx k Vfi [x] x r ~ 1 (35) 

with x n = 1 — r n = 0.427 and xk = 1 — rx = 0.954. The functions Z U \T] are the (r — l)-momenta of 
the photon-neutrino kernels and can be calculated analytically by using Eqs. ([31~]) and ([33]) . They are 
shown in Fig. [3] for where the dotted, dashed and solid lines corresponds to v = v^, V^ and (f M + 17^), 
respectively. The red line is the approximate expression Z v +„ \T] = 0.71 — 0.16 (r + 0.1) obtained 
by Eqs. (5,6) of jTHj where the relationship between 7-ray and neutrino spectra produced by simply 
parametrized primary proton spectra was studied. One sees that our calculation predicts ~ 10% larger 
neutrino fluxes. This is not surprising considering that the calculation of [18] (which is based on [29] ) 
do not include kaons (and assume equal production of tt + , ir~ and 7T°)0 

The photon spectrum is expected to have a cutoff at high energies. Thus, we can write <3? 7 = 
N £ 7 [E] E , where N is a normalization factor and the adimensional function £ 7 [E] which modulates 
the photon spectrum is normalized to 1 at low energies. If 7-rays originate from hadronic processes, 
the cutoff cannot be too sharp, since the spectral features of the parent CR spectrum are diluted in 
hadronic cascades, as was pointed out in [22]. If the function £ 7 [i£] is sufficiently smooth, we can extract 
it from the integrals in Eqs. ([30]) and ([32]) . obtaining: 



$ V [E] =NZ u [T]UE]E~ r (36) 



with 



&„ [E] - ^r^T ( °- 380 x - r ^7 [E/x K ] + 0.013 x K r ^ [E/x K ] + £ 7 [E/x m ] [ dx ^ [x] x 
^i/ M [J-J V Jo 

Cf m [E] ~ — ^r ( 0.278 x n r C 7 [E/x n ] + 0.009 x K r ^ [E/x K ] + ^ 7 [E/x m ] [ dx k^ [x] x 

L ~VuV-\ \ JO 



(37) 



where x m = 0.59 approximately coincides with the maximum of the function k v [x] x r_1 when T ~ 2 — 3. 
The above relations shows that, beside being suppressed by a factor Z U [T] with respect to the photon 
flux, the neutrino fluxes are also shifted to lower energies. For T = 2 the numerical factors of the various 
terms in the r.h.s. of Eqs. (j37p become 0.069, 0.012, 0.130 and 0.051, 0.008, 0.137 respectively, showing 
that the dominant contribution is provided by neutrino produced in muon decays (i.e., last terms in 
the r.h.s. of Eqs. ([3?]) ). As a consequence, the features in the photon spectrum, such as the presence of 
a cutoff at an energy EU C , are essentially reproduced in the neutrino spectrum at an energy lower by 
a factor E VtC /Ej tC = x m ~ 0.59, as also noticed by [18]. 



6 We do not compare with the results of the recent publication [34] because they do not provide sufficient details to 
reproduce their results. In particular, we are unable to reproduce the line labeled "Vissani, 2006" in their Fig. 2 which, 
in the intention of the authors of [34], should amount to an application of |16| . 
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2.4 

r 



Figure 3: The junction Z U \T] for v^ (dotted), v^ (dashed) and v^ + v^ (solid). The red line is the approximate 
expression Z v +■$ [r] = 0.71 — 0.16 (r + 0.1) obtained by Eqs. (5,6) of U8f . 



5.2 Using raw data 

RXJ1713. 7-3946 is presently the best studied SNR. It has been observed by H.E.S.S. during three 
years from 2003 to 2005 [Sj. The data expands over three decades, exploring the energy interval 
Ey = 0.3 — 300 TeV. The energy resolution of the experiment is equal to about 20% and the photon 
spectrum is sampled in 25 bins SE^/E^ ~ 0.2 plus three larger bins at high energy. In a previous 
paper [22], we have discussed the information that the observational data provide on the parent CR 
spectrum in the SNR. Here, we use the observational data to calculate the v^ and V^ fluxes emitted by 
this object and to estimate the event rate expected in neutrino telescopes. The proposed procedure, 
which does not require any parametrization of the photon flux, is general and can be applied to any 
other 7-transparent source. 

As a first step, we "rescale" the photon and neutrino fluxes according to: 



Vl [E] = 
ME] = 



$y[E] 
$»[E] 



E a 
E a 



(38) 
(39) 



where v = z/„, v^. For a proper choice of the parameter a, the "rescaled" fluxes are expected to vary 
slowly with energy. In the following, we adopt: 



a 



2.5 



(40) 



which is particularly appropriate for RX J1713. 7-3946 (see, e.g., Fig. 3 of [22] and related discussion). 
We indicate with ipj ± Aipj the (rescaled) photon flux measured in the j—th. energy bin, centered 
at a photon energy Ej and covering the energy range (Ejj n f,Ej tSup ). We can approximate the photon 
flux by: 

^[E 7 ]=J2<fjW j [E y ] (41) 

3 

where Wj-[£L] are rectangular functions which describes the various energy bins (i.e., Wj[Ey] = 1 for 



E 



jM 



< £L < E. 



J, sup 



and zero elsewhere). We immediately obtain from Eqs. (|30I32|) the relation: 



<Pv„ [E] 



J2vjWj[E] 
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0.5 1 1.5 

Log[E v /l TeV] 



Figure 4: The v^ (solid) andV^ (dotted) flux expected from the SNR RX J 1713. 7-3946 according to the H.E.S.S. 
ry—ray data. The shaded area are obtained by propagating the observational uncertainties in the "f-ray data. The 
red dotted lines correspond to the atmospheric neutrino flux in the vertical direction integrated over an angular 
window 0.5° (lower line) and 1.0° (upper line). 



^i E ] = ^2<PjWj[E\ 



(42) 



where: 



Wj [E] = 0.0453 Wj [E/(l-r 1t )]+ 0.0116 Wj [E/(l - r K )] + / — k u \x] x a WAE / x] 

Jo x 

/•l dx 
Wj [E] = 0.0331 Wj [E/(l - tv)] + 0.0080 Wj [E/(l - r K )] + —k Vfi [x] x a Wj [E/x] (43) 

jo x 

Rels. (|42|) give the neutrino fluxes as a linear combination of the observational values tpj of the 
photon flux. The functions tUjf-E 1 ] and u7i[-E] describe the contribution that each data point gives to the 
reconstructed neutrino flux at the energy E. The uncertainty in the neutrino fluxes can be evaluated 
by propagating linearly the observational errors Aipj. We obtain: 






(44) 



Similarly, the correlation between the values of the u^ flux at two different energies can be calculated 
by: 

Y jk ^lw k [E]w k [E'] 



Qv, 



\E,E'} 



Ej^jWjlE^JEi^fw^E] 



(45) 



Analogous expressions can be obtained for z/„ flux by replacing ifj[-E] — ► Wj[E] in the above expressions. 

The results of the proposed method are displayed in Fig. |4] where we show with solid (dashed) line 

the ffj, flux (v^ flux) expected from RX J1713. 7-3946 according to the H.E.S.S. observational 7— ray 

data. The H.E.S.S. data cover the energy range £L = 0.3 — 300 TeV. In this interval we have described 



the photon flux according to Eq. (|41j) and we have not relied on theoretical assumptions. It is, however, 
unnatural to assume that the photon spectrum vanishes outside the probed region. For this reason, 
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we have continued the photon spectrum at low energy (EL < 0.3 TeV) assuming $ 7 = I (E /ITeV)"™ 
where f3 = 2 and / = 2.47 • 10~ 11 /(cm 2 s TeV). The assumed low-energy photon spectrum smoothly 
connects with the low energy bins observed by H.E.S.S experiment (see, e.g., [8], [22])- In order to make 
minimal assumptions we have not extrapolated the photon spectrum at high energies. The relevance 
of the assumed low and high energy behavior can be studied by considering different extrapolations for 
the gamma ray flux outside the region probed by H.E.S.S. experiment. In this way, we have been able 
to verify the spectral region directly constrained by the data is E v ~ 0.3 — 100 TeV. 

We can see that the neutrino and antineutrino spectra are well described by a power law with 
spectral index 7 ~ 2 at low energies and a cutoff/transition region at E v ~ 3 — 5 TeV. The ratio 
3>j7 /$!,„ is nearly constant and equal to about 0.93. The shaded areas describe the observational errors 
in the neutrino fluxes and are obtained by propagating the uncertainties in the gamma ray data. The 
neutrino fluxes are well constrained at low energies where the relative uncertainty is at few per cents 
level. The information degrades at high energy. We have A<& u /& u ~ 30% at E v = 30 TeV and much 
worse at larger energies. We remind, for completeness, that one has to consider also a systematic error 
in our calculation equal to about 20% due to uncertainties in hadronic cross section, neutrino oscillation 
parameters and the approximations implicit in our method, as conservatively estimated in the previous 
sectionHI 

The red dotted lines in Fig.[4]show the atmospheric neutrino flux which provides a diffuse background 
for SNR neutrino detection. The relevance of this background is, clearly, reduced if one is able to observe 
the source with a good pointing accuracy. The lower (upper) red dotted line is obtained by integrating 
the muon neutrino flux in the vertical direction given by [36] over an a angular window 0.5° (1.0°), 
which correspond to the angular resolution of neutrino telescopes at energies equal to about 1 TeV. We 
remind that, below E u ~ ITeV, the angular response of the detector is mainly determined by kinematic 
of the detection process and, thus, cannot be improved. We see that the spectral region in which the 
signal is expected to be larger than the background is above ~ ITeV quite close to cutoff in the v 
spectrum. The red lines in Fig. H] should be intended as lower limits of the atmospheric v background. 
The atmospheric neutrino flux, in fact, depends on the zenith angle and increase by about one order 
of magnitude in the horizontal direction. Moreover, the angular dimension of the galactic SNRs can be 
comparable with the detector angular resolution. The SNR RX J1713. 7-3946, as an example, subtends 
an angle in the sky equal to about 1.0°. 

In the next section, we will calculate the event rate produced by RX J1713. 7-3946 in a telescope 
located in the Mediterranean sea. We conclude this section by discussing whether we can expect 
neutrino fluxes much larger than that from RX J1713. 7-3946. In this respect, it was noted in [22] that 
the 7-ray spectrum of RX J1713. 7-3946, corresponds to a flux of CR protons with a energy equal to 
about 0.05 x 10 51 erg (consistent with Ginzburg & Syrovatskii hypothesis) colliding with a molecular 
cloud of approximately 300 solar masses, which represents a quite favorable situation. In other words, 
it difficult to imagine CR fluxes in SNRs much more energetic than this or a much larger target mass. 
SNRs can be closer than RX J1713. 7-3946 which has a distance from us equal to about D ~ lkpc. 
Keeping the same parameters, this would increase the 7 and v fluxes as D~ 2 and, correspondingly, 
the expected signal in v telescopes. However, there will be no gain with respect to the atmospheric 
v background, if the linear dimension L of the sources is the same as RX J1713. 7-3946, since the 



7 The wiggly behavior of the neutrino spectra is not physically significant. It reflects the statistical fluctuations of the 
photon data in the simple interpolation scheme proposed in Eq. (|41[) . To reduce the effect, we have applied a Gaussian 
smearing to the neutrino fluxes predicted by Eq. (|42p on scales smaller than 8E V / E v ~ 0.03. This is equivalent to smooth 
the photon flux on scales SE-y/E-, < 0.03 and does not erase any significant features. The smoothing scale is, in fact, 
much smaller than the H.E.S.S. energy resolution (equal to about 20%) and than the size of the energy bins. 
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Figure 5: Left Panel: The function n Vll [E u ,E t } 1 }E~ a (solid) and n T: ^[E l/ ,E t \ i ]E 1 a (dotted) for a = 2.5 
and E th = 0.05, 0.2, 1, 5, 20 TeV. Right Panel: The function n[E n , E th \E- a for a = 2.5 and E th = 
0.05, 0.2, 1, 5, 20 TeV. 

angular dimension 6 = L/D increase with decreasing distances. For this reason, the best candidates 
are close and young SNRs which had no time to expand to large radii and which are also expected to 
produce spectra which extend to larger energies, having the most energetic proton less time to escape. 
A promising candidate could be the young and close SNR named Vela Jr (angular size 2°), as suggested 
by the extrapolations of the first H.E.S.S. observations [T7 1 l22 | 121] . 

6 Events in neutrino telescopes 

Charged-current interactions of v^ and V^ produce muons and antimuons that can be observed in neu- 
trino telescopes. The number of muons and antimuons reaching an area A and in a time of observation 
Tis: 

/■oo 

AW = hvA-T- dE ^ [E] x Y u [E, E th ](l - a^ [E] ) + (^ - F M ) (46) 

where E is the energy of the neutrino at the point of interaction and E t ^ is the energy threshold for 
muon detection. Following [14]: (1) we assume A = 1km and T = 1 solar year; (2) we consider that 
the SNR RX J1713. 7-3946 is visible for /u v = 78% of the sideral year, as in Antares location; (3) the 
neutrino absorption coefficient a v (E), averaged over the daily location of the source, is calculated for 
standard rock; (4) the probability Y^E^E^) that a neutrino of energy E produce a muon that reach 
the detector with an energy larger than _E t h is calculated for water. 
The above equation can be rewritten in the form: 



N, 



H+fi 



-E-th 



dEu 
EL 



n 



,A\E v ,Eto\$ Vli [E v } + (v» 



(47) 



where the function: 



n„ M [E v , E th ] = f hv -A-T- Y U [E U , E th ] ■ (1 - a^ [E v ]) ■ E v 



(48) 



gives the relative contribution of neutrinos with energy E v to the total event rate above a detection 
threshold E t ^. In the left panel of Fig. [5] we show the function n u [E u , E t \^\E~ 2 ^ as a function of E u for 
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the selected values Eth = 0.05, 0.2, 1, 5, 20 TeV. One sees that the bulk of the signal is expected to 

arise from the spectral region E v ~ 0.3 — 100 TeV which is directly constrained by photon observational 

data, provided that the detection threshold is below few TeV. 

The previous point can be much better appreciated by rewriting the expected event rate in v- 

telescopes in terms of the observed photon flux. By using the results of the previous sections, we 

obtain: 

f°° dE 
N n+-P= ^n[£ 7 ,£ th ]$ 7 [£ 7 ] (49) 

JE th E 1 

where the function n[E.y,Eth] is given by: 

r E ~i dE 
n[E^E th \ = / — ^ K^EjE^n^E^E^] + (i/„ - F M ) (50) 

and i<T° sc [x] is the photon-neutrino kernel given in Eq. (|19p . The function n[E^, E t h]E~ 2 ' 5 is shown as 
a function of E~ in the right panel of Fig. [SJ We see immediately that the spectral region probed by 
H.E.S.S. data (i.e., E 7 = 0.3 — 300 TeV) is the most relevant energy region to derive expectations for 
future neutrino telescopes when the detection threshold is E& < 1 TeV. In this case, neutrinos that 
effectively contributes to the v— telescopes signals are connected to photons which are observed with 
Cerenkov 7— ray telescopes. If the threshold is much larger than this, on the contrary, the function 
n[Ej, E t h]E~ 2,5 is peaked outside the probed region and one is forced to rely on theoretical extrapola- 
tions. 

By following a procedure analogous to that adopted in the previous section, we can calculate the 
expected signal from RX J1713. 7-3946 directly from photon observational data. If the photon flux is 
approximated as in Eq. (|4"Tj) . we obtain: 

N H+-p = J2 i Pj Jl j[ E th] (51) 

3 

where <pj = 3>j • E^ is the observational value for the (rescaled) photon flux in the j—th. energy bin and 

the factors: 

/dP 
-g 2 n[£ 7 , E th \ E~ a Wj [E y ] (52) 

"weights" the contribution of each energy bin to the signal above the threshold E t ^ (Wj[E^] = 1 for 
-E^inf < E,y < Ej jSup and zero elsewhere). The uncertainty in the predicted neutrino signal can be 
evaluated by propagating linearly the observational errors A(fj, obtaining: 



AAW= /;L A ^^th] 2 (53) 

The number of neutrino events from SNR RX J1713. 7-3946 in Antares location after one year of 
data taking is reported in Tab. [2j We see that a km 2 class neutrino telescope will be able to collect 
few events per year, if the threshold will be lower than about 1 TeV. The observational uncertainty is 
less than 10%, indicating that H.E.S.S. data very well constrains the expected neutrino signal. The 
above estimates do not include possible contributions from the low (_E 7 < 0.3 TeV) and high energy 
(E~ > 300TeV) tails of the photon flux, which are not constrained by observations. However, these 
contributions are expected to be negligible (at the few per cents level), if the threshold is lower than 
few TeV. We remind that our calculation is affected by ~ 20% systematic error due to uncertaintie in 
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Table 2: The number of neutrino events, N^ + -p, from SNR RX J 1713. 7-3946 expected in Antares location per 
km 2 per year for various energy thresholds, according to H.E.S.S j-ray data. The uncertainty of the predicted 
signal, AN^+j;, is obtained by propagating H.E.S.S. observational errors. The atmospheric neutrino background, 
is estimated from \36\j (see text for details). 



yvrAtmo 



E th (TeV) 


%+H 


AAW 




jyAtmo 


0.05 


5.65 


0.35 


0.06 


20.5 


0.2 


4.67 


0.33 


0.07 


6.6 


1 


2.44 


0.28 


0.11 


1.1 


5 


0.57 


0.17 


0.30 


0.1 


20 


0.08 


0.07 


0.95 


0.007 



hadronic cross sections, neutrino oscillation parameters and the approximations implicit in our method, 
as conservatively estimated in the previous sections. 

The above results show that the detection threshold will be a crucial parameter. Clearly, it cannot 
be much larger than 1 TeV, otherwise the event rate would be negligible. However, it cannot be much 
lower than this, due the presence of the atmospheric neutrino background. At present, the atmospheric 
neutrino fluxes have a quite large uncertainty due to imprecise knowledge of primary CR spectra at 
earth, of hadronic cross sections and the poor understanding of charmed particles production (see [36]). 
However, in the next future, they will be measured by neutrino telescopes. 

In order have feeling of the background levels, we approximate the atmospheric neutrino fluxes 
according to $>£ tmo [E v ] = <& Atmo [E^] = C ■ (E^/lTeV) -3 - 6 , as it is appropriate in the energy range of 
our interest. The normalization constant C = 1.8 x 10~ n TeV~ cm" 2 s" 1 has been obtained by fitting 
the vertical neutrino flux of [36] (which underestimates the average atmospheric neutrino flux in the 
direction of the SNR) and integrating over an angular window = 1°. By using this parametrization, 
we obtain the number of events given in the last column of Tab. [2j Clearly, these values are purely 
indicative. They allow us, however, to conclude that there is no real gain to lower the threshold below 
the TeV level, in agreement with the finding of [37]. In conclusion, the best choice seems E t h ~ 1 TeV, 
since it allows us to probe neutrino emission: i) in an energy region very-well constrained by 7-ray 
data; ii) collecting few events per year; Hi) with a signal-to-background ratio order one or larger. 



7 Summary 

Future observations in neutrino telescopes have the potential to test an important theoretical paradigm, 
that the young SNRs are the main site of acceleration of the galactic cosmic rays. At present, it is 
useful to use the existing observations of very high energy gamma rays from SNRs to derive quantitative 
predictions for neutrinos from these sources. The main results obtained in this paper are the following: 
i) We have discussed a conceptually and computationally simple method to extract precise pre- 
dictions for neutrinos from supernova remnants and from other hypothetical sources, transparent to 
their gamma rays. This method (that is based and that elaborates on our previous work on the sub- 
ject [16l HH [221 [21]) is superior to other ones present in the literature, in that it does not need a 
preliminary parametrization of the gamma ray observations and, as we demonstrated, permits one to 
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use directly the gamma ray data as an input. 

ii) Our method allows us to propagate easily the 7-ray observational errors and, thus, to understand 
how well neutrino fluxes and the signal in future neutrino telescopes can be constrained by 7-ray 
observations. Our analysis, and specifically the application to the best studied SNR (RX J1713.7- 
3946) shows that the present, successful program of observations of very high energy gamma rays from 
certain supernova remnants (in particular, with imaging arrays of Cerenkov telescopes) covers the right 
energy region to derive expectations for the forthcoming neutrino telescopes. As an example, the signal 
produced in ^-telescopes by RX J1713. 7-3946 is predicted with an observational uncertainty equal to 
about ~ 10%. 

Hi) We have discussed the sources of theoretical and systematic uncertainties in our calculation. In 
particular, we have checked that the error arising from unknown neutrino oscillations parameters is very 
small (at the level of ~ 2%), as a result of the partial cancellation of the (much larger) anti-correlated 
contributions of electron and muon neutrino oscillation probabilities to the total error budget |f| The 
main source of uncertainty in the predictions is due to the modeling of the hadronic interactions and was 
conservatively estimated at ~ 20% level. This error could be more precisely assessed after a systematic 
comparison of existing numerical codes and overview of the available experimental data. 

iv) It seems possible, at least for the best observed SNR (RX J1713. 7-3946), to succeed and detect 
a neutrino signal with exposures of the order of yearx km 2 , provided that the detection threshold in 
future neutrino telescopes will be equal to about ~ ITeV. Another promising and perhaps better young 
SNR is Vela Jr, whose higher part of gamma ray spectrum (above 20 TeV) is still to be studied. Due to 
the presence of the atmospheric neutrino background, it does not seem really much useful to lower the 
threshold for neutrino observation below the TeV region. This is in agreement with the finding of [37j . 

In summary, we have shown that it is possible to obtain precise and reliable predictions of the 
neutrino signals expected in neutrino telescopes from SNRs. The observation of neutrinos from these 
objects would amount to a proof of the existence of the expected over-density of CR and, thus, to a 
confirmation of the hypothesis that young SNRs are the main site of acceleration of galactic CR. 

Acknowledgment 

This work was partially supported by the High Energy Astrophysics Studies contract number ASI- 
INAF 1/088/06/0, by the MIUR grant for the Projects of National Interest PRIN 2006 "Astroparticle 
Physics" and by European FP6 Network "UniverseNet" MRTN-CT-2006-035863. We thank F. Aharo- 
nian, P. Blasi, M.L. Costantini, P.L. Ghia, P. Lipari, F. Lucarelli and G. Riccobene for useful discussions. 



3 This is also due the vicinity of a physical boundary, almost saturated by oscillations, as can be seen from in Fig. [2] 
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